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Abstract 

Humanoid robots locomote by making and breaking contacts with their 
environment. A crucial problem is therefore to find precise criteria for a 
given contact to remain stable or to break. For rigid surface contacts, the 
most general criterion is the Contact Wrench Condition (CWC). To check 
whether a motion satisfies the CWC, existing approaches take into account a 
large number of individual contact forces (for instance, one at each vertex of 
the support polygon), which is computationally costly and prevents the use 
of efficient inverse-dynamics methods. Here we argue that the CWC can be 
explicitly computed without reference to individual contact forces, and give 
closed-form formulae in the case of rectangular surfaces - which is of practi¬ 
cal importance. It turns out that these formulae simply and naturally express 
three conditions: (i) Coulomb friction on the resultant force, (ii) ZMP in¬ 
side the support area, and (iii) bounds on the yaw torque. Conditions (i) and 
(ii) are already known, but condition (iii) is, to the best of our knowledge, 
novel. It is also of particular interest for biped locomotion, where undesired 
foot yaw rotations are a known issue. We also show that our formulae yield 
simpler and faster computations than existing approaches for humanoid mo¬ 
tions in single support, and demonstrate their consistency in the OpenHRP 
simulator. 


1 Introduction 

From the viewpoint of the robot, establishing contact with the environment amounts 
to constraining a certain number of Degrees Of Freedom (DOF) of an end-effector 
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link. For instance, a sliding planar contact imposes three equality constraints (two 
on the orientation of the link and one for the link-to-surface distance) while a fixed 
contact constraints all six DOFs of the end-effector link. A stability criterion can 
be seen as a set of inequality constraints describing the conditions under which 
these equalities are preserved. 

The venerable Zero-Moment Point (ZMP) [Q criterion may be the best-known 
example of stability criterion for rigid surface contact. It is known 13 to be nec¬ 
essary but not sufficient for contact stability: in particular, it says nothing about 
possible sliding in the horizontal plane and rotation around the vertical axis (yaw). 
Yet, humanoid robots are often subject to significant yaw moments during single 
support phases, resulting in undesired foot rotations. Attesting the importance of 
this problem, recent works have been using upper-body motions to compensate for 
these yaws while continuing to use ZMP EldJ. 

A more principled (and general) way to address this problem is to consider in¬ 
dividual contact forces distributed on the contact surface, as can be found e.g ., in 
bipedal balance control 0 and motion planning (6] 0 1. This approach yields a 
stronger stability criterion than ZMP, and accounts for both the sliding and yaw 
rotation. It is however hampered by redundancy : the vector of contact forces has 
many more components (three times the number of contact points) than the degree 
of the contact constraint (six). This redundancy makes the resolution of the equa¬ 
tions of motion fundamentally harder, as illustrated by the fact that state-of-the-art 
Inverse Dynamics based on QR-decomposition 0 only apply to non-redundant 
variables^ In Time-Optimal Path Parameterization (TOPP), redundancy prompted 
the use of further contact approximations J9J or expensive polytope projection al¬ 
gorithms [0. In the present paper, we argue that such workarounds are no longer 
necessary if one uses the correct contact representation (for instance, both 0 and 
101 boil down to a single matrix inversion for a biped in single-support, as we will 
see in Section [5]). 

The key insight here is that the condition that individual contact forces lie in 
their respective friction cones can be replaced by the condition that the contact 
wrench belongs to a certain wrench cone ifTOl . The contact wrench naturally solves 
the redundancy issue, as its dimension is minimal (six). It was advocated as a gen¬ 
eralization of ZMP in 0, along with a stability theorem, and applied to walking 
pattern generation on rough terrains [11, H3- However, this theorem makes the 
same “sufficient friction” assumption as ZMP, which means the resulting criterion 
does not account for sliding and yaw rotations. Besides, the contact wrench is 

1 When this is not the case, the authors advocate the use of Singular Value Decomposition to 
compute new independent variables; however, it is unclear how to compute the inequality constraints 
applying to these new variables. 
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computed from individual contact forces, which we argue yields unnecessarily ex¬ 
tensive calculations: following iFTOi the wrench cone can be computed explicitly 
from the sole geometry of the contact surface. 

In this paper, we derive the closed-form formulae of the wrench cone in the case 
of rectangular support areas, which is of practical importance since most humanoid 
robot feet can be adequately approximated by rectangles. This result helps simplify 
dynamics computations, as we will see for humanoid motions. It also provides 
an analytical description of “yaw friction”, from which we derive a simple and 
principled control law to prevent undesirable yaw rotations. 

The rest of the paper is organized as follows. In Section |2| we recall the defi¬ 
nitions related to contact stability. In Section [3} we discuss the physics of surface 
contact and give a theoretical justification for the practice of considering individual 
contact forces at the vertices of the contact polygon. Then, in Section |4| we derive 
a closed-form expression of the wrench cone in the case of rectangular contact ar¬ 
eas. We apply the resulting solution in a humanoid experiment in Section [5]bef ore 
concluding in Section [6j 


2 Background 

2.1 Contact Forces and Contact Wrench 

Consider a robot with n degrees of freedom making N point contacts with the 
environment, at points Ci,..., Cn in the laboratory reference frame (RF). The 
equations of motion of the robot are: 

N 

M(q)q + h(q,q) = S T r a + ^J(C' i ) T f i , (1) 

2=1 

where q, q, q are the n-dimensional vectors of DOF values, velocities and accelera¬ 
tions, M is the n x n inertia matrix, h(q, q) the n-dimensional vector of gravity and 
Coriolis forces. In case the robot has n a actuated joints, r a is the n a -dimensional 
vector of torques at the actuated joints and S is a n a x n joint selection matrix. 
Finally, for each i E [1, N], fj is a 3-dimensional vector of contact force and J (Cf) 
is the 3 x n translation Jacobian calculated at point C x . 

We assume that both the environment and the contacting link are rigid bodies. 
Thus, interactions between them can be represented by a single contact wrench 
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W = (f, r), which we can compute from contact forces as: 


f d = Y'fi, ( 2 ) 

i 



i 


where O is the origin of the link RF. The contact jacobian J wr is the 6 x n matrix 
obtained by stacking vertically J (O), the translation Jacobian computed at O, and 
J ro t, the rotation Jacobian of the link, both taken with respect to the absolute RF. 
With these definitions, we have the following property: 

N 

4rW = ]Tj(a) T fi. (4) 

i —1 

(See Appendix |A| for a proof.) Then, the equations of motion can be rewritten as: 

M(q)q + h(q,q) = S T r a + 4 r W. (5) 


2.2 Contact Stability 

Assume now that the robot is in a given state (q, q). The accelerations q and 
generalized forces exerted on the robot (actuated torques or contact forces) are 
bound by a complementarity conditions [13]. For the sake of the explanation, let 
us consider first the simple case of a single translation coordinate x, as depicted in 
Figure [I] 


X 



Figure 1: Block on a horizontal surface with one DOF. 

Under Coulomb’s friction model, either of the following situations occurs: 

• Fixed contact: x = 0 and the contact force obeys \fl\ < pf™, where f l c and 
are respectively the horizontal and vertical components of f c and p is the 
static friction coefficient; 
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• Sliding contact: x > 0 and the contact force obeys f l c — —pkfc> where 
is the kinetic friction coefficient. 

The acceleration x and contact force f c are in a complementary relationship: when 
one is equality-constrained, the other is inequality-constrained. Similarly, for the 
general case of a 6 -DOF end-effector in contact, the translational and rotational ac¬ 
celerations of the link are in a complementary relationship with some generalized 
contact forces 7 (here, 7 = W or (fi,..., f at)). The contact mode describes, for 
each variable in a complementary relationship, whether it is equality- or inequality- 
constrained. In contact stability, we are interested in th t fixed contact mode where 
the position and orientation of the end-effector are equality-constrained to a refer¬ 
ence value. 

Definition 1 (Weak Contact Stability) A contact is weakly stable when there ex¬ 
ists a solution (q, r a , 7 ) of the equations of motion satisfying the fixed contact 
mode for all contacting links. 

That is to say, 

• for each contact (z), the relative velocity and acceleration at contact are zero: 
Jwrq = 0 and Jwrq = - J^q, 

• actuated torques r a are within torque limits, 

• complementary forces 7 satisfy their inequality constraints (friction cones 
or wrench cone). 

This formulation has been widely used in the literature. In approaches based on 
inverse dynamics, the conditions on (q, q) are first enforced kinematically, then 
torques and complementary forces are solved 0 [ 8 ]|. 

The “weakness” in the definition above refers to the notions of strong and weak 
stability, as stated by lfl3ll . Strong stability happens when all solutions to the equa¬ 
tions of motion satisfy the fixed contact mode. Note that choosing between contact 
forces and the contact wrench changes the equations of motion (Equations ([I]) and 
([5]), respectively), but the underlying contact stability is the same by Equation ([4]). 
In the rest of the paper, we will always refer to contact stability in the weak sense. 

3 Surface Contact 

Suppose we take contact forces fi,..., f/v as complementary variables to the po¬ 
sition and orientation of the contacting link. Let ff and denote the normal and 
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the tangential components of the contact force f \. Coulomb friction provides the 
complementary inequalities: 


• Unilaterality: 

fi> o, 

(6) 

• Non-slippage: 


(7) 


where y is the static coefficient of friction. 

However, when the contact is done through a surface and not through a set 
of points, the reality of contact is that of continuum mechanics. In this case, the 
action of the environment at the contact surface S is described by two quantities: 
a scalar field p(x,y) corresponding to normal pressure , and a two-dimensional 
vector field cr(x,y) for tangential mechanical stress. Figure [2]-( A) illustrates these 
two fields for a rectangular contact area. For convenience, we also will denote by 
v = cr(x, y) + p(x, y) n, where n is the unit vector normal to the contact surface. 
The equations of motion become 


Mq + h = S T r a + J J(C xy ) T v(x,y)dxdy (8) 

where J (C xy ) is the 3 x n translation Jacobian calculated at the point of coordinate 
C xy on the surface (taken in the laboratory RF). The wrench resulting from v(x,y) 
is 


f = f J^v(x,y)dxdy, (9) 

r = f J Ocjy A v{x, y)dxdy, (10) 


where O is the origin of the link RF. Under Coulomb friction, the inequality con¬ 
straints for v(x,y) are: 


• Unilaterality: 

• Non-slippage: 


p(x,y) > 0, 
<r(x,y) || < np(x,y). 


( 11 ) 

( 12 ) 


Note that taking a constant // in Equation ( fl2l ) is an approximation: in reality, the 
friction coefficient n(x,y) varies with the position on the surface. 

In the present literature, surface contact is often modeled using sets of contact 
points: sometimes more than required by positional constraints ( e.g ., nearly 30 per 
contacting link in 0) or one at each vertex of the convex hull (e.g., in 0). The 
proposition below gives a theoretical justification for the latter practice. 
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Proposition 1 Assume that the contact surface S is a convex polygon with ver¬ 
tices Ci,..., Cat- If there exists i/(x, y)( x ^ y )^s satisfying complementary inequal¬ 
ity constraints Cy}-©’ then there will exist contact forces applied at ... ,Cjst 
and summing up to the same contact wrench that satisfy complementary inequality 
constraints @-Q- 



(A) 



Figure 2: Contact in the surface plane. (A) Example of stress/pressure fields. Red 
discs indicate the magnitude of pressure by their size. Blue arrows show tangential 
stress. (B) Notations used in Section[4] 

Proof: consider pressure and stress fields summing up to W. By convexity, one 
can find strictly positive functions ai(x, y),..., ctk( x i y) such that ai{x , y) = 
1 and each point C xy G S can be written C xy = y)Ci. Then, define for 

each vertex C % a force 

fj := oti{x, y)v(x, y)dxdy, 

By positivity of the afs, it is straightforward to check that all ff > 0 and ||f/|| < 
pff. In addition, this expression of f \ ensures that the resulting wrenches are equal, 
i.e. 9 0 = © and @ = ([TDJ. □ 

We furthermore argue that, when the friction coefficient /i is assumed to be 
constant as in Equation ( fl2| ), the two conditions are equivalent. The complete 
proof of this equivalence, which involves reconstructing pressure and stress fields 
given local and boundary conditions, is however out of the scope of this paper. The 
bottom line of this argument is that using contact forces at vertices of the convex 
hull completely describes the dynamics of surface contact. The wrench cone that 
we derive in the next section will share the same characteristic. 
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4 Wrench Cone for Rectangular Surfaces 

Consider a rectangular area (C1C2C3C4) as depicted in Figure [ 2 ] (B). We calcu¬ 
late the contact wrench W at the center O in the link’s RfJ^] Let us denote by 
(/f, ff , f z ) the three components of the contact force at point C\ Under the com¬ 
mon linear approximation of friction cones, Coulomb inequalities become 


1 /fU/fl < 

fi > 

The wrench cone is then given by the following proposition. 


0 


(13) 

(14) 


Proposition 2 There exists a solution (fi,..., £4) satisfying inequalities (fl 3 ])-(p 4 |) 
if and only if there exists a wrench W = (f x , f y , f z ,r x ,r y ,r z ) such that: 


in 

< 

vf z 

(15) 

\p\ 

< 


(16) 

r 

> 

0 

(17) 

\t x \ 

< 

Yf z 

(18) 

\l~ y \ 

< 

Xf 

(19) 

7min — 

T Z 

fi Tnax 

(20) 


where 


T min = ~V(X + Y)f Z + \Yf X - ^T X \ + \Xfy - VT V \ 
W d = +v(X + Y)f z -\Yf x +fiT x \-\Xfy + vT y \ 


The proof of this proposition is given in Appendix |A.1| The validity of this 
expression was also tested empirically with a script available at fl4ll . 

Let us now detail each line of the wrench cone. The first two inequalities <d- 
( [T6| ) correspond to the usual Coulomb friction. Inequalities ( fTT] ), ( fl8] ) and ( fl9| ) are 
equivalent to the ZMP condition. The last inequality ( [20] ) provides a bound on the 
admissible yaw torque that was implicitely encoded in the contact-force model. 
Note how this relation is more complex than a mere “no rotation occurs while r z 
is small enough”, as it is coupled with all other components of the contact wrench. 
Notably, the “safest” value is not zero but: 


def / 

Aafe vi 


min “t - 7"max)/2 

= sgn(-/ 3: r a: )min(y|/ a: |, y u|r a: |) 

+ sgn(-f y T y )mm{X\f y \,/j,\T y \), 


2 multiply by the link’s rotation matrix for a wrench in the absolute RF 
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Figure 3: Snapshots of the retimed motion. The motion lasts 7.285 s. Time stamps 
are shown below each frame. This motion stresses all components of the wrench 
cone. The first part stresses the pitch by moving the COM forward. The second 
segment extends the arms, then keeps the left arm still while moving the right arm, 
thus stressing the roll. At the same time, the chest pitch is actuated back and forth, 
which stresses the yaw. Finally, the waist performs an elliptic motion (back and 
forth, up and down) throughout the whole motion, thus stressing the translation of 
the contact foot. 

where sgn is the sign function. From ( [20] ), t z may deviate from r sa f e by at most 

n(X + Y)f - max(Y\f x \, h\t x \) - max(X\f y \, h\t v \). 

We see that higher tangential forces or roll-pitch torques reduce the range of ad¬ 
missible yaw torques. In particular, when these other constraints are saturated (for 
instance when the ZMP reaches a corner of the support polygon), T sa f e is the only 
solution that prevents the contact from breaking. Therefore, r z = r sa f e appears as 
a sensible control law to prevent undesired yaw rotations. 

5 Experiment 

We test the validity of the Contact Wrench Condition (CWC) in the integrated 
simulator OpenHRP with a model of the HRP4 robot. Note that OpenHRP has its 
own contact model where forces are distributed along the edges (not corners) of 
the contact surface. 

We implemented the CWC condition within the Time-Optimal Path Parame¬ 
terization framework (TOPP), a well-known projection of system dynamics along 
a pre-defined path that has been used for motion planning of humanoid robots 
01130. We considered the case of a single contact at the left foot and designed 
a motion that would challenge all six contact DOFs. In single contact, the contact 
wrench is fully determined by the unactuated rows of the equation of motion, i.e., 

W = (PJwr)- 1 ? [M(q)q + h(q, q)] 
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Figure 4: Zoom on the contact forces computed by OpenHRP at the left foot. The 
first row corresponds to the retimed motion (7.29 s), the second to the 10% uni¬ 
form acceleration of the retimed motion (6.55 s) and the third to the 15% uniform 
acceleration (6.19 s). Time stamps are shown below each column. The retimed mo¬ 
tion maintains surface contact with contact forces distributed all around the contact 
area, except between 4.2 and 4.8 s where they are mostly on the left half of the foot. 
The 10%-accelerated version goes to point contact (3.0 s) and line edge contact (3.6 
to 4.8 s) but does not fall. At 15% acceleration, the CWC violation overcomes the 
stabilizer’s ability and the robot falls. 

where P = I — S is the unactuated line selection matrix. Note how computing the 
constraint projector is straightforward when using the contact wrench: it is not the 
case when contact forces are used, as we observed in previous research j9l (where 
we resolved the force redundancy with a force-binding model). 

For this experiment, we designed by hand a set of eleven key postures. The 
geometric path was obtained by interpolating Bezier curves between these postures 
and applying a kinematic filter to fix the position and orientation of the support foot 
on the ground. To get a feasible trajectory from this path (i.e., to compute the timing 
information) we used the open-source TOPP solver fl6!l . The solver takes as input 
the path and a vector representation of the system dynamics along it, which is easy 
to compute once one knows the projectors mentioned above (see (9) for details). 
Figure [3] shows a timelapse of the final retimed motion. Videos are also available 
online at fi4l . 

An interesting thing to note about TOPP is that, because of time optimality, the 
retimed motion always saturates at least one of the contact constraints. In an ideal 
setting with perfect system dynamics, the obtained trajectory should therefore be 
at the limit of contact: it would execute correctly, but the contact would break as 
soon as one tries to perform the motion faster at any time instant. We observed this 
phenomenon in the experiment, as depicted in Figure [4] 

There are, however, a few points to discuss. In order to ensure proper tracking 
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of the computed trajectory, we used HRP4’s stabilizer module, which occasionally 
changes the actual joint angles. In turn, this behavior violates the assumption made 
by TOPP that the robot follows a precise geometric path, and the constraints may 
get under- or over-saturated. To alleviate for this issue, we added safety margins 
in TOPP’s conditions: we scaled the contact area by 45% and set a smaller friction 
coefficient /i = 0.4 (versus fi = 0.8 in OpenHRP). Consequently, it was possible 
in practice to accelerate the retimed motion by about 5% reducing the total 
duration by 5% with uniform timescaling) and obtain a successful execution. How¬ 
ever, we still observed the expected phenomenon with relatively small changes: for 
a 10% acceleration (by uniform timescaling), the humanoid started to lose surface 
contact, but the stabilizer was still able to recover; for a 15% acceleration, the 
violation of the contact constraint was too large and the robot fell. See Figure [4] 


6 Conclusion 

In this paper, we calculated the closed-form expression of the wrench cone for 
rectangular contact surfaces. This formula has several implications. First, it is very 
simple, making computations much easier than any previous formulation based 
on contact forces. Second, it describes concisely the phenomenon of yaw friction 
by a double-inequality that is, to the best of our knowledge, novel. From these 
bounds, we derived a simple control law to avoid undesired foot rotations, a recur¬ 
ring problem for bipeds in single contact. Finally, we showed how our criterion can 
give simpler and faster computations than contact forces for humanoid motions in 
single support, and demonstrated it with dynamic motions simulated in OpenHRP. 
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A Proof of Equation [4] 


The rotation Jacobian satisfies u = J ro tq, where u> is the rotation velocity of the 
link. An interesting consequence of this property is that for any point C on the link 
and any vector u, one has 



( 21 ) 


Next, the position of any point C L of the link in the absolute RF is related to that 


of its origin O by C{ — O + OCi. The corresponding translation Jacobian is 
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The first term of the expression equals the translation component of J wl . W. By 
applying ( |2lT ), we see that the second term equals J rnt (()(A f, ). Factoring J mt 
out of the summation yields the rotation component of J wr W. □ 


A.l Calculation of the Wrench Cone 

The wrench is defined by ([T|)-([T|i as: 

f X = ff + P + P + ff 
P = fi + fl + /| + fl 
P = fi + /I + /I + /I 

r x = Y(ff - /f - /| + /|) 

TV = -X(f{ + /I-/I-/I) 

tZ = X(ff + /|-/|-/|)-F(/f-/|-/| + /*). 

By unilaterality ( fl4l ) we have f~ > 0, so we can define: 


Ki 

■= P/pP 

Ci 

:= r x /Yf z 

k 2 

■= P/pP 

c 2 

: = T y/xf z 

k 3 

:= t z /h{X + Y)P 

A 

■■= P/EiP 

Px 

:= X/(X + Y) 

Py 

:= Y/{X + Y) 

a x 

L l 

■■= P/pP 

CP? 

1 

■■= ff/pff 


and normalize the system by dividing each row accordingly. From the non-slippage 
constraint we have af , a v { and Di G [—1,1]. Then, introduce the new vari¬ 
ables: 

lx = Oi\D x + Oi\D± = ol\ D 2 + of £>3 

7?/ = cii-Di + at y 2 D 2 7y = a|£>3 + 

We can reduce the complete system in two ways. First, using the fact that the rela¬ 
tion M from a to 7 is a linear surjection from [—1, l ] 8 to M[—1, l ] 8 = {\^ %| < 
Dj + £&} (the computation of antecedents being straightforward). Since there is 
no other constraint on the c^’s than their domain and relation to 7 ^’s, one can ob¬ 
tain an equivalent system by replacing a G [—1, l ] 8 by 7 G M[—1, l] 8 . Then, the 
three equations in £^’s are: 

1 = £1 + £2 £ £3 £ £ 4 ? 

Ci — D\ — £2 — £3 + £> 4 , 

C2 = —D\ — £2 + -£3 + £4? 

By linear combination, we can use them to rewrite M[— 1, l ] 8 as: 2\^/ x \ < 1 + C\, 
2\i x \ < 1 — Ci, 2 | 7 y | < 1 — C 2 , 2 \^' y \ < 1 + C 2 . Finally, using the same equations, 
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one can express all D{ s as functions of, e.g ., D 4 . The inequalities D G [0, l ] 4 
become: 


-1 + 

Cl 

< 

2D 4 

< 

l + C’i 

C\ + 

c 2 

< 

2D 4 

< 

2 + C\ + C 2 

-1 + 

C 2 

< 

2D 4 

< 

l + c 2 

0 


< 

2D 4 

< 

2 


This system has solutions if and only if all lower bounds are smaller than all upper 
bounds. Matching all pairs of lower and upper bounds one by one (we will show an 
example of this technique below in a more complex situation), one can check that 
all these inequalities boil down to C\ G [—1,1] and C 2 G [—1,1]. The complete 
system is now: 

Kl = 1x + V x K 2 = 7y + Vy 

K 3 = Pxily - iy) - Pyilx - l' x ) 

2|7*l < l + C'i 2\i x \ < 1-Ci 

2173/1 < i-c 2 2 | 7 ;i < 1 + C 2 

And (Ci, C -2 ) € [—1, l] 2 . One can use the first three equations to eliminate the 
redundant variables 7 ^, 7 ^ and 7 ^, expressing them as functions of 7 ^ in the in¬ 
equality constraints. After simplification, the resulting system is: 


^Pylfx 

< 

Pyi 1 + 

Cl) 


( 22 ) 

^Py^x 

< 

Py(l - 

Cl) + 2p y K\ 


(23) 

2 Pylfx 

< 

Px( 1 — 

C 2 )-K 3 +p y K 1 - 

PxK 2 

(24) 

2 Pylx 

< 

Px( 1 + 

C 2 ) — K 3 + PyKi + 

PxK 2 

(25) 

^Pylx 

> 

-Pj/l 1 

+ Ci) 


(26) 

ZPylx 

> 

“Pl/U 

— Cl) + 2pyK\ 


(27) 

ZPylx 

> 

-ja(i 

-C 2 )-K 3 +p y K 4 

~P X K 2 

(28) 

^Pylx 

> 

-ja(i 

+ C 2 )-K 3 +p y K 4 

+ p x k 2 

(29) 


And (Ci,C 2 ) G [—1, l] 2 . There exist a solution 7 ^ if and only if all of its lower 
bounds are smaller than all of its upper bounds. Let us match all pairs of lower 
boundsi 


and upper bounds f22])-([25]). One can check that: 

• ([26j < {22} Ci > -1 

• {26} < {23} Kx > -1 


• {26]) < {24]) K 3 - p y Ki + p x K 2 - p y C 1 + p x C 2 < 1 

• {26} < {25} K 3 - p y Ki - p x K 2 - pyCi - p x C 2 < 1 
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• {27} < {22} 4=> Ki < 1 

• {27} < {23} 4=> Ci < 1 

• {27) < {24) K 3 + p y Ki + p x K 2 + p y Ci + p x C 2 < 1 

• {27}<{25}^ K 3 + PyKl - p x K 2 + PyCl - p x C 2 < 1 

• ([28j < ([22J -K 3 + p y Kx - p x K 2 - p y Ci + p x C 2 < 1 

• ([28} < ([23} —K 3 - pyKx - p x K 2 + p y Ci + p x C 2 < 1 

• ([28} < ([24} C 2 < 1 

• ([28} < ([25} A 2 > -1 

• |29} < ([22} -K 3 + p y Ki + p x K 2 - pyCi - p x C2 < 1 

• ([29} < ([23} -X 3 - Pyifi + p x K 2 + pyCi - p x C 2 < 1 

• ([29} < ([24} i^ 2 < 1 

• ([29} < ([25} C 2 > -1 
Consequently, the complete system becomes: 


k 3 

< 

1 - PyKi - 

PxK 2 - PyCi - 

PxC 2 

k 3 

< 

1 - PyKi + 

PxK 2 - PyCi + 

PxC 2 

k 3 

< 

1 + PyKi - 

p x K 2 + PyCi - 

PxC 2 

k 3 

< 

1 + PyKi + 

PxK 2 + PyCJ \ + 

PxC 2 

k 3 

> 

-1 + PyKi 

+ P X K 2 - PyCi 

~PxC 2 

k 3 

> 

-1 + PyKi 

- PxK 2 - PyCi 

+ PxC 2 

k 3 

> 

-1 - PyKi 

+ PxK 2 + PyCi 

~PxC 2 

k 3 

> 

-1 - PyKi 

- PxK 2 + PyCi 

+ PxC 2 


And (i^i, K 2 , Ci, C 2 ) E [—1, l] 4 . In a more concise form, these last eight inequal¬ 
ities can be written K% ^ —1 H-pyl-fC — Ci| p x \K 2 — C 2 | and K% ^ + 1 -Pyi-^i + 
Ci | p,r l\ > ■ C> . We conclude by de-normalizing all inequalities. □ 
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